Magnetic phase separation in EuB 6 detected by muon spin rotation 
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We report results of the first muon-spin rotation measurements performed on the low carrier 
density ferromagnet EuBg. The ferromagnetic state is reached via two magnetic transitions at 
T m = 15.5 K and T c = 12.6 K. Two distinct components are resolved in the muon data, one 
oscillatory and one non-oscillatory, which arise from different types of magnetic environment, and 
we have followed the temperature dependence of these components in detail. These results provide 
evidence for magnetic phase separation and can be interpreted in terms of the gradual coalescing of 
magnetic polarons. 

PACS numbers: 76.75.+i, 75.47.Gk, 75.50.Cc 
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Europium hexaboride has attracted recent interest be- 
cause it exhibits colossal magnctoresistance (CMR)* and 
it has been suggested that its semiconductor-semimetal 
transition results from the overlap of magnetic polarons 2 . 
EuBg crystallises into a simple cubic structure (space- 
group PmSm) with divalent Eu ions ( 8 <SV/2) at the corners 
of the unit cell and Bg-octahcdra at the body-centred po- 
sitions, and is a ferromagnet at low temperatures^. Spe- 
cific heat and magnetization measurements reveal that 
this state is reached via two distinct transitions at T m = 
15.5 K and T c = 12.6 K 2 ^. Neutron diffraction measure- 
ments show that a small spontaneous magnetic moment 
begins to grow on cooling below T m , but does not become 
significant until T c is reached, below which point the 
moment shows a more usual mean-field like behaviour—. 
The magnetic ordering is accompanied by a sharp drop 
in the resistivity which is strongly field dependent and 
gives rise to a large negative magnetoresistance 2 . This 
transition from a semiconductor at high temperatures 
to a semimetalSi&S (or possibly a self-doped compen- 
sated semiconductor*^) at low temperatures is reminis- 
cent of the metal-insulator transition seen in manganitcs 
exhibiting CMRi*. Detailed measurements of resistivity 
and magnetization^ show that this transition is associ- 
ated with T m . It is thought that magnetic polarons could 
be responsible for this behaviour, with the upper mag- 
netic transition and drop in resistivity caused when the 
bound carriers overlap and percolate, and the lower tran- 
sition caused by a true transition to a bulk ferromagnetic 
state 2 ". Further support for this explanation comes from 
the observation of polaronic features, possibly associated 
with itinerant holes*^, below ~ 30 K in Raman-scattering 
spectral. The trapping of carriers to form bound mag- 
netic polarons provides an explanation for the upturn 
in resistivity observed on cooling through 30 K 2 *^ and 
the negative magnetoresistanc o 12 ' 14 . However, a direct 
observation of overlapping polarons in EuBg has so far 
been elusive. 

In this paper we present the results of /iSR experi- 
ments on EuBg which not only provide further evidence 
for the two distinct magnetic transitions but are also able 



to resolve two components arising from muons stopping 
in two different types of environment below T m . The 
first component is an oscillating signal and can be asso- 
ciated with muons that stop in a locally ferromagnetic 
environment. The second component is a Gaussian sig- 
nal and arises from the muons that come to rest in a 
paramagnetic environment. This provides clear evidence 
for magnetic phase separation below T m . 

Our /iSR experiments were carried out using the 
DOLLY instrument at the Paul Scherrer Institute (PSI) 
in Switzerland and the DEVA beamline at the ISIS pulsed 
muon facility in the UK. In our fiSK experiments, spin 
polarised positive muons (/i + , mean lifetime 2.2 ^s, mo- 
mentum 28 MeV/c) were implanted into polycrystallinc 
EuBg. The muons stop quickly (in < 10~ 9 s), without 
significant loss of spin-polarisation. The time evolution 
of the muon spin polarisation can be detected by count- 
ing emitted decay positrons forward (f) and backward 
(b) of the initial muon spin direction due to the asym- 
metric nature of the muon decayi*. In our experiments, 
positrons are detected using scintillation counters placed 
in front of and behind the sample. We record the num- 
ber of positrons detected by forward (JVf) and backward 
(iVb) counters as a function of time and calculate the 
asymmetry function, G z (t), using 



G z {t) 



N t (t) - a cxp N h (t) 



N f (t) 



(1) 



where cv CX p is an experimental calibration constant and 
differs from unity due to non-uniform detector efficiency. 
The quantity G z (t) is then proportional to the average 
spin polarisation, P z (t), of muons stopping within the 
sample. The muon spin precesses around a local mag- 
netic field, B (with a frequency v = (^ fi /2i:)\B\, where 
7^/2^ = 135.5 MHz T" 1 ). 

Examples of asymmetry spectra measured at PSI are 
shown in Fig.^ There are three distinct temperature re- 
gions. In the lowest temperature data (T < 10 K), there 
are clear oscillations in the measured asymmetry, demon- 
strating that the sample does indeed make a transition 
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FIG. 1: Muon decay asymmetry plots for EuBg at different 
temperatures. The solid lines are fits of the data to Eq. 
The data for T = 11.5, 14 and 15 K are fitted with a non-zero 
Gaussian amplitude, whilst the others are not. The data were 
taken at PSI. 



to a locally magnetically ordered state. This oscillating 
signal is superposed on a slow exponential relaxation un- 
observable in Fig. ^ but visible at longer times. In the 
second region (10 K < T < T m ), the amplitude of the os- 
cillatory component decreases and a Gaussian component 
appears, with a decay rate that decreases as the temper- 
ature is increased. In the third region (T > T m ), the 
amplitude associated with the Gaussian component de- 
creases until, when above 16 K, only the slow exponential 
relaxation remains. The need for the fast relaxing Gaus- 
sian term in the fits is clearly motivated by the topmost 
panel of Fig. ^ where the purely exponential relaxation 
seen at 17 K is compared with the results for a tempera- 
ture of 15 K, which is just inside the intermediate region. 

In order to best follow these changes, the data over 
the whole studied temperature range were fitted to the 



function 

G z (t) = A exp exp + A>s C exp (-jr) cos(27rz/t) 

+ A gauss exp(-cr 2 i 2 ) + A hgl (2) 

where Ab g represents a time-independent background 
due to muons stopping in the silver that surrounds the 
sample, Ti and T2 are the longitudinal and transverse re- 
laxation times, and A cxp , A osc and A gauss are the ampli- 
tudes of the exponential, oscillating and Gaussian terms 
respectively. Plots of the fitted parameters against tem- 
perature are shown in Fig. [21 Vertical lines drawn at 
T m and T c approximately divide the plots into the three 
temperature regions discussed above. In Fig. ^e) the 
asymmetry amplitudes of the three different terms are 
shown, which provide the main motivation for this divi- 
sion. The exponentially relaxing fraction is seen to drop 
sharply from the maximum value at T m ; this is expected 
in a polycrystalline sample when a transition into an or- 
dered state occurs. In the lowest temperature region the 
remaining amplitude is accounted for by the oscillating 
fraction, so that here the sample is uniformly magneti- 
cally ordered, but at intermediate temperatures there is 
also a contribution from the fast relaxing Gaussian frac- 
tion. 

The fact that the asymmetry amplitude is shared be- 
tween the two fast signals in the intermediate region is 
evidence that the muons are stopping in two different 
kinds of environment, that is, magnetic phase separation 
is occurringiSiit The amplitudes of the two signals are 
expected to be proportional to the volume fraction of 
each phase. Phase separation is important between T m 
and T c , but we find it disappears below T c , in contrast 
to the predictions of Ref. 

A Gaussian relaxation can result from a field at the 
muon site which is static but randomly distributed in 
magnitude. The parameter a is related to the width of 
the field distribution as a 2 = ^y tl 2 (B 2 )/2. Fig.^d) shows 
that a increases as the temperature is lowered in the 
second region, indicating that the field distribution at 
the muon sites is becoming wider as the sample becomes 
more ordered at lower temperatures. 

Fig-I3 a ) shows that the oscillations develop below T m , 
and their frequency v increases fairly smoothly as the 
sample is cooled through T c rising to a maximum of ~ 
67 MHz at low temperature (corresponding to a field at 
the muon site of ~ 0.5 T). The frequency v is propor- 
tional to the magnetization and the data in Fig. Ufa) 
match well with the temperature dependence of the Eu 
magnetic moment measured with neutron scattering 3 . 
Only one muon precession frequency is observed, strongly 
suggesting that there is a single set of equivalent muon 
sites in the structure. There are two candidate sites. 
The first is at the centre of a Bg~ octahedron, and the 
second is at the face-centres of the unit cell (in the cen- 
tre of the shortest B-B bond which is between atoms 
in adjacent unit cells). Using the Eu moment measured 
previously 3 and assuming a low temperature magnetic 
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FIG. 2: Temperature dependence of the parameters deter- 
mined from fits of Eq. |5] to EuBe asymmetry spectra. The 
panels correspond to (a) the oscillation frequency, u, (b) and 
(c) the relaxation rates 1/Ti and I/T2, (d) the Gaussian frac- 
tion width parameter, a, and (e) the amplitudes for each of 
the three components used in the fits. The dashed lines in (e) 
are guides for the eye. The vertical dashed lines indicate the 
positions of T c and T m . 



structure with the moments pointing along the [111] di- 
rection, the dipole field (-Bdip) can be calculated at both 
these possible sites. At the centre of a boron octahe- 
dron, the dipole field cancels by symmetry. The face- 
centre positions are all magnetically equivalent and yield 
7/J-Bdip|/27r = 144 MHz. Additional contributions to 
the field at the muon site arise from the Lorentz field 
(/i M sat /3 = 0.42 T, corresponding to 57 MHz), the de- 
magnetization field and the hyperfine contact field, and 
preclude a definitive assignment of the site. 

The longitudinal relaxation rate, 1/Ti, reflects the dy- 
namics of the fields being probed. For rapid fluctuations, 
1/Ti oc llY< q \&B{q)\ 2 T{q), where \SB(q)\ is the ampli- 
tude of the fluctuating local field and r(q) is the Eu-ion 
correlation time at wavevector q. In the paramagnetic 
phase the spin fluctuations are so rapid that the mea- 
sured relaxation rate is small. As the sample is cooled 
and the critical region is approached, the correlation time 
becomes longer and the relaxation rate rises and peaks 
close to T m (Fig. Gib)). In contrast, I/T2 (which is pro- 
portional to the width of the ordered field distribution 
corresponding to the oscillating fraction) rises dramat- 
ically on cooling through T c and subsequently falls on 
further cooling (Fig. Etc)). This emphasises that the or- 
dered and fluctuating fractions in the sample are distinct 
and follow different temperature dependences. 

These observations fit in well with the polaron per- 
colation picture for the intermediate temperature range; 
muons stopping in ferromagnetic regions of overlapping 
polarons give rise to the oscillating signal, while most 
muons stop in the intervening paramagnetic regions. In 
these paramagnetic regions the Eu-ion correlation time 
is short, so the muons are not depolarized by the local 
fluctuating Eu moments, but by the distribution of fields 
that result from the nearby ferromagnetic regions. As the 
sample becomes more ordered the paramagnetic regions 
shrink, with a corresponding effect on the field distribu- 
tion. 

In order to study the slower relaxation in the high tem- 
perature region in more detail, data were collected at the 
ISIS facility. An example is shown in Fig.Qfa) and clearly 
shows the presence of both a slow and a fast relaxation 
rate. Therefore the data were fitted to the function^ 



G z (t) = A f exp(-Aft) + A s exp(-A s i) + A t 



(3) 



where Af tS and Af. s are the amplitudes and relaxation 
rates of the fast and slow exponential components (note 
that Af represents relaxation due to fast dynamics and 
hence slow spin relaxation). The non- relaxing back- 
ground fraction A^ g was held fixed. Af and A s were 
found to be approximately constant with temperature, 
with Af about half the value of A s , and were fixed in the 
fits. The relaxation rates themselves appear temperature 
independent above the transition, and it was possible to 
keep Af fixed at 0.08 MHz (see Fig. |3J). On cooling to- 
wards T c , A s begins to decrease, until only one relaxation 
rate can be resolved below 9K. This is possibly because 
the two rates become very similar in magnitude, and the 
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FIG. 3: Panel (a) shows a typical asymmetry spectrum mea- 
sured at ISIS. The solid line is a fit to Eq. [2] and the dot- 
ted lines show the form of the two exponential terms used. 
Panel (b) shows the two relaxation rates measured at ISIS, 
Af , s (filled symbols) , and the longitudinal relaxation rate mea- 
sured at PSI, 1/Ti (open symbols). The solid line is a fit of a 
T 2 ln(T/A) function to the PSI rate for temperatures below 
14 K. 

presence of Af may mask any further drop of A s . 

The relaxation rate A s matches well with the tem- 
perature dependence of 1/Ti, as shown in Fig. Efb), 
and the two can be identified with each other. At low 



temperatures, a contribution to 1/Ti can be fitted by 
1/Ti oc T 2 In (T/A) (see Fig.^b)), appropriate for scat- 
tering by two-magnon processes in ferromagnetsi?. How- 
ever, the component due to fast fluctuations produces too 
slow a relaxation to be convincingly included in the fits 
to the PSI data. Nevertheless, the observation of two 
relaxation rates above T m allows us to infer the presence 
of spatial inhomogeneity (a very similar effect has been 
found in Lao.67Cao.33Mn03 with /iSRi£), which persists 
even at 115 K. Above T m , the temperature dependence 
of Af and A s is weak and featureless. Any polaron dynam- 
ics are presumably too fast to be followed by the muon. 
The dramatic changes observed below T m can therefore 
be attributed to a large change in the time-scale of po- 
laron dynamics, such as might be expected when polarons 
overlap. Although the polaron volume fraction changes 
little at a percolation transition, their arrangement could 
be changed so that the average size of a single polaron is 
larger and its dynamics are much slower. 

In conclusion, /iSR measurements have allowed us to 
follow the very unusual development of ferromagnetism 
in EuBg through the transitions at T m and T c from a 
local viewpoint. These results reveal two distinct and 
spatially separate regions of the material associated with 
different magnetic behaviour. Such magnetic phase sepa- 
ration is qualitatively consistent with a picture based on 
coalescing polarons. 
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